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We use scanning tunneling spectroscopy to investigate the filled and empty electronic states of 
superconducting single-unit-cell FeSe deposited on SrTiO3(001). We map the momentum-space 
band structure by combining quasiparticle interference imaging with decay length spectroscopy. 

In addition to quantifying the filled-state bands, we discover a F-centered electron pocket 75 meV 
above the Fermi energy. Our density functional theory calculations show the orbital nature of empty 
states at F and explain how the Se height is a key tuning parameter of their energies, with broad 
implications for electronic properties. 


The extraordinary potential of interface engineering to 
generate novel electronic properties is exemplified by a 
single unit cell (lUC) of FeSe deposited on SrTiOs [1], 
which exhibits an order-of-magnitude increase in its su¬ 
perconducting transition temperature (Tc up to 110 K 
[2]) compared to bulk FeSe (Tc = 9.4 K [3]). Not only 
does this finding elevate the Tc of iron-based supercon¬ 
ductors (Fe-SCs) above the liquid nitrogen temperature, 
it also opens the door to designing Fe-SC/oxide het¬ 
erostructures with novel phases and yet higher Tc. A 
key to understanding and realizing these phases is a com¬ 
plete measurement of the electronic structure of filled and 
empty states. 

Electronic band structure is pivotal in determining the 
pairing symmetry of Ee-SCs. The generic Eermi surface 
(ES) of Ee-SCs consists of electron pockets at the Bril- 
louin zone (BZ) corner M and hole pockets at the zone 
center F [4] . A prevalent spin-fiuctnation model suggests 
that repulsive antiferromagnetic excitations of wave vec¬ 
tor (tt, tt) can give rise to pairing between the electron 
and hole pockets if the order parameter reverses sign, re¬ 
sulting in S-\ _superconductivity [5, 6]. However, in lUC 

FeSe/SrTiOs, the F hole pocket sinks entirely below the 
Fermi energy (Ep) due to electron doping [3]. This chal¬ 
lenges the 5^_picture; nevertheless, functional renor¬ 

malization group (FRG) calculations have shown that 
electronic bands lying within the spin fluctuation energy 
scale below Ep can still influence the pairing channel. In 
fact, the energy of the sunken F hole pocket is predicted 
to toggle the relative stability between sign-preserving 
S++ and sign-changing d pairing symmetries [8, 9]. 

A natural question is whether low-lying bands above 
Ep can similarly renormalize the effective interaction. In 
general, the landscape of empty states in Fe-SCs remains 
largely unexplored by experiment. A full band struc¬ 
ture mapping is particularly crucial in lUC FeSe/SrTiOs, 
where in addition to the usual Coulomb repulsion and 
spin fluctuations, even higher energy phonon modes may 


be at play [9-11], and the magnitudes of their energy 
scales relative to the neai-Ep bands determine the su¬ 
perconducting ground state. 

Here we map the multiband electronic structure of 
lUC FeSe/SrTiOs by two complimentary scanning tun¬ 
neling microscopy (STM) techniques: (1) quasiparticle 
interference (QPI) imaging [13] and (2) decay length 
spectroscopy [14]. In the first technique, impurity scat¬ 
tering of quasiparticles generates interference patterns 
with characteristic dispersive wave vectors q{uj) that can 
be inverted to reconstruct the band structure. Since q 
is the momentum transfer, QPI imaging resolves only 
relative momentum coordinates between two states. In 
the second technique, the absolute, in-plane momentum 
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FIG. 1. (color online) (a) Typical topography of in-situ-grown 
FeSe/SrTiOs. Setpoint: 4 V, 5 pA. (b) Line cut along the ar¬ 
row in (a). The inset illustrates the underlying crystal struc¬ 
ture. (c) Atomically-resolved topography of single-unit-cell 
(lUC) FeSe/SrTiOa. Setpoint: 50 mV, 250 pA. (d) dl/dV 
spectrum of lUC FeSe/SrTiOs, T = 4.3 K. Bias oscillation 
Frms = 0.7 mV. 
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FIG. 2. (color online) Quasiparticle interference imaging, real space, (a) Topography (setpoint: 50 mV, 500 pA) and (b-e) 
conductance maps g{Y^uj) (setpoint: 100 MQ, V-ms = 1-4 mV) of a 20 nm x 20 nm field of view with as-grown defects. Images 
were drift-corrected following Ref. [1]. 


of quasiparticles can be extracted from the decay of 
their tunneling current with increasing sample-tip sepa¬ 
ration. By combining the two momentum-resolved tech¬ 
niques, we discover a F electron pocket 75 meV above 
Ef- Our density functional theory (DFT) calculations 
reproduce the presence of empty states at F, and fur¬ 
thermore explain how their energies are tuned by the Se 
height hse- 

We grew films of FeSe on Nb-doped SrTiO3(001) 
(0.5%) via molecular beam epitaxy (MBE). The sub¬ 
strates were pretreated with deionized water for 90 min 
at 80°C, followed by an 02-anneal for 3 h at 1000°C. 
We then transferred the substrates into our MBE cham¬ 
ber (base pressure 1x10“^^ Torr) and degassed them 
at 670° C. We deposited EeSe by co-evaporating Ee 
(99.995%) and Se (99.999%) with a molar flux ratio of 1:6 
and substrate temperature 520°C. Afterwards, we typi¬ 
cally annealed the samples for an additional 2 h between 
500-600° C before transferring them through ultra-high 
vacuum to a home-built STM for imaging at ^4.3 K. 

Eigure 1(a) shows a typical film topography, with re¬ 
gions of bare SrTiOs and lUC or 2UC of EeSe. We dis¬ 
criminate these regions based on their terrace heights. 
Erom the line cut in Eig. 1(b), we observe a 3UC SrTiOs 
step to be 1.19 ± 0.05 nm (bulk c-axis lattice constant is 
0.3905 nm [15]), the 1-2UC EeSe step to be 0.57 ± 0.05 
nm, and the bare SrTiOs-lUC EeSe step to be 0.34T0.02 
nm (all measured at 4 V sample-tip bias). We will 
hereafter focus on the lUC EeSe terraces. Eigure 1(c) 
presents an atomically-resolved topography of lUC EeSe, 
with lattice constant a = 3.9 A. Each bright spot corre¬ 
sponds to a surface Se atom in a Se-Ee-Se triple layer. A 
representative dljdV spectrum on a clean area exhibits 
a gap of A = 14 meV (Eig. l(ci)), similar in magnitude 
to other reports of superconducting gaps in this material 
[3, 16]. We note appreciable spectral inhomogeneity in 
lUC EeSe/SrTiOs, but further study is needed to quan¬ 
tify its correlation with substrate disorder. 

To image QPI, we acquired conductance maps 
g{r^uj) = dl/dV{r^eV) over fiat regions of lUC 
EeSe/SrTiOs with moderate concentrations of as-grown 
defects (Eig. 2(a)). Several energy maps of one repre¬ 


sentative region are presented in Eigs. 2 (b-e), display¬ 
ing clearly dispersive interference patterns. To identify 
the momentum-space origin of the scattered quasipar¬ 
ticles, we compared the Eourier transform amplitudes 
\g{q,uj)\ to simulated autocorrelations of the spectral 
function A{k^uj) = [17]. Eor sim¬ 

plicity, we used the bare Green’s function G~^{k,uj) = 
uj + iS — £a{k)^ with parabolic bands £a{k) and broaden¬ 
ing (5 = 5 meV. The main result is presented in Eigs. 3(a- 
i), which compare \g{q^uj)\ to theoretical predictions for 
three representative energies. We discuss each in turn: 

uj = 10 meV, Eigs. 3(b,e,h): Close to we observe 9 
ring-like intensities in \g{q^uj)\^ centered about reciprocal 
lattice vectors G = (0, 0), (±27r/a, 0), (0, ±27r/a), and 
(±27r/a, ±27r/a). These intensities arise from scattering, 
modulo G, within electron Eermi pockets at the zone 
corner M (labeled 1 in Eig. 3) [3]. 

UJ = —66 meV, Eigs. 3(c,f,i): Sufficiently below Ep^ 
we observe additional scattering channels pointing to the 
emergence of the F hole pocket seen by angle-resolved 
photoemission spectroscopy (ARPES) [3]. Intrapocket 
scattering between F pockets is labeled 2 in Fig. 3, while 
interpocket scattering between F and M pockets is la¬ 
beled 1-2 in Fig. 3. 

cj = 80 meV, Figs. 3(a,d,g): Above Ep^ we discover a 
third pocket. Intrapocket scattering (labeled 3 in Fig. 3) 
is clearly resolved in \g{q^uj)\^ but interpocket scatter¬ 
ing with the M electron pockets (expected intensity at 
{'Kja^Tila) modulo G) appears to be suppressed. In gen¬ 
eral, the autocorrelation of A{k^uj) yields the set of all 
possible scattering channels, but more complex theories 
that encode spin [18] or orbital [19] selectivity in the scat¬ 
tering T-matrix are needed to explain their relative in¬ 
tensities. In this case, the empirical suppression of F-M 
scattering leaves some ambiguity as to the absolute mo¬ 
mentum {k) location of the new pocket. 

To visualize the full QPI evolution, Fig. 3(j) shows an 
azimuthally-averaged intensity plot of |^(g^,ci;)|, where 
Qr is measured relative to G = (27r/a, 0) as shown in 
Fig. 3(h). In total, we observe three dispersing branches: 
two electron-like (labeled 1 and 3) and one hole-like (la¬ 
beled 2). Branches 1 and 2 correspond to a M electron 
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FIG. 3. (color online) Quasiparticle interference imaging, momentum transfer (q) space, (a-f) Theoretical simulations, A(k,(ju) 
and its autocorrelation, for three representative energies, (g-i) Fourier transform amplitudes \g{q,uj)\ of conductance maps 
(four-fold symmetrized for increased signal), (j) Azimuthally-averaged intensity plot of \g{qr,uj)\, where Qr is measured relative 
to G = (27r/a,0). The superconducting gap is marked by 2A. 


pocket and a F hole pocket, while branch 3 awaits fur¬ 
ther identification. A parabolic fit to branch 1 over the 
given energy range in Fig. 3(j) yields an effective mass 
enhancement m*/m = 2.0T0.1 and a carrier concentra¬ 
tion of 0.08 e~ per Fe from a Luttinger count, assuming 
a degenerate pocket [3, 20]. 



FIG. 4. (color online). (a,b) Energy dependent decay length 
A(ct;), extracted from exponential fits to the tunneling current 
as the tip is retracted from the sample at a fixed bias (inset 
schematic). Fits were performed in the current range [10 pA, 
500 pA], two of which (250 meV and 50 meV) are shown in 
(a). Dashed horizontal lines indicate calculated values of A at 
the F and M. (c) dl/dV spectrum. Erms = 2.8 mV. Vertical 
lines mark extrema in the numerical derivatives of A(ct;) and 
dlldV. 


To determine the absolute momentum k of QPI branch 
3, a complimentary momentum-resolved STM technique 
is needed. Here we utilize decay length spectroscopy 
[4, 14, 21], a general tool which allows the full recon¬ 
struction of /c-space band structure from STM. Tersoff 
and Hamman [23] showed that a sample state of in-plane 
momentum k\\ has density which decays towards the vac¬ 
uum with length A given by 

(2A)2 ^ ’ 

where ^ is the average of the sample and tip work func¬ 
tions. Figures 4(a,b) show the energy dependent decay 
length A(co’), extracted from exponential fits to the tun¬ 
neling current as the sample-tip distance is increased at 
a fixed bias. Near the sample states have large mo¬ 
mentum near M and smaller decay length. Below Ep: 
a steep increase in A(ci;) accompanies the onset of a hole 
pocket at F, as states with low momentum become avail¬ 
able for tunneling. The fact that a similar rise in A(ci;) 
occurs above Ep indicates that branch 3 in Fig. 3(j) is 
also located at F. If we interpret the large-jcjl value of 
A = 0.462 ±0.001 A as arising from states with k ^ 0 ^ we 
find ^ = 4.46 ± 0.03 eV from Eq. (1), then we can com¬ 
pute the expected A(co’) = 0.318 ± 0.001 A for energies 
where the only states come from momenta near M. In¬ 
deed, the measured A(co’) at small \uj\ closely matches the 
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FIG. 5. (color online) (a) Band structure of free-standing 
single-unit-cell (lUC) FeSe, calculated in the generalized gra¬ 
dient approximation (GGA). Structural parameters: lattice 
constant a = 3.90 A, Se height hse = 1-45 A. (b) Energies of 
the five F bands shown in (a) vs. hse (the band represented by 
orange is degenerate). The GGA values of hse for lUG FeSe 
(a = 3.90 A fixed) and bulk FeSe (a = 3.68 A relaxed) are 
marked, as well as the experimental value for bulk FeSe. The 
Fermi energy Ep expected from electron doping is marked in 
(a,b). (c) Gharge density isosurfaces (yellow) at = 0 for 
the five F bands, shown in two perspectives. The histograms 
depict the orbital compositions. 

expected value of X{\k\ = \f2'Kjo). Step-like features as¬ 
sociated with the onsets of these pockets are also detected 
with dl/dV spectroscopy (Fig. 4(c)). From extrema in 
the numerical derivative d?I/dV^^ which closely match 
those of dXjduj (vertical shaded guides in Figs. 4(b,c)), 
the band edges of the F hole and electron pockets are -65 
meV and 75 meV. 

A consistent band structure for lUC FeSe/SrTiOs is 
now established, comprising M electron pockets span¬ 
ning Ep and F hole and electron pockets lying below 
and above Ep. For further insight, we use DFT to 
compute the band structure of free-standing lUC FeSe 
via the generalized gradient approximation (GGA) [24] 
and projector augmented wave (PAW) method as im¬ 
plemented in the Vienna Ab-Initio Simulation Package 
(VASP) [25, 26]. We use a BZ sampling of 9x9x1 and 
an energy cutoff of 450 eV. We apply Methfessel-Paxton 
smearing [27] with a = 0.1 eV. Figure 5(a) shows the cal¬ 
culated bands with structural parameters a = 3.90 A, hse 
= 1.45 A. Due to electron doping, Ep should be adjusted 
to intersect only the M pockets. Typical band renormal¬ 
ization factors range from 4-5 in lUC FeSe/SrTiOs [28], 
but for the qualitative discussion that follows, we do not 
rescale the bands. 


Experimentally, hse is unknown. Simulations show 
that the binding geometry of lUC FeSe/SrTiOs varies 
with Ti 02 oxygen deficiency, which creates electropos¬ 
itive sites that distort Se positions [29]. Without mi¬ 
croscopic knowledge of the buried interface, we calculate 
band structures for a range of hse values and track the 
energies of the F bands (Fig. 5(b)). While all bands shift 
slightly, the lowest-lying F electron pocket in Fig. 5(a) 
undergoes a pronounced monotonic decrease in energy 
with increasing hse- Figure 5(c) shows the charge den¬ 
sity isosurfaces at k = 0 and orbital compositions for each 
band. Only the lowest-lying F electron pocket carries sig¬ 
nificant Se 4p character in addition to Fe 3d character, 
so it is most affected by the Fe-Se distances. The charge 
density plot suggests an antibonding configuration of Fe 
3 dx 2 _y 2 and Se Apz orbitals, which explains the increase 
in pocket energy with greater overlap of Fe and Se states. 
Our calculation reveals a crucial connection between hse 
and empty electronic states. 

Previous reports have predicted that Se/Te heights 
tune the Fe exchange constants in iron chalcogenides and 
hence the magnetic order [30], which is oddly absent in 
FeSe [31] and unknown in lUC FeSe/SrTiOs. Here, we 
discuss another implication of hse- As seen in Fig. 5(b), 
the F electron and hole pockets cross at large values of 
hse- Recently, Wu et al. have proposed that nontrivial 
Z 2 topology may be realized in lUC FeTei_a,Sea, [32]. In 
particular, when the gap between the F electron and 
hole pockets falls below 80 meV, spin-orbit coupling can 
invert the bands. We measure to be 140 meV from 
Fig. 4(c); thus, lUC FeSe/SrTiOs could possibly lie in 
proximity to a topological phase transition. 

In summary, we have quantified both the filled and 
empty state band structure of lUC FeSe/SrTiOs, and 
discovered a new F-centered pocket emerging around 75 
meV above the Fermi level. Our work has several impor¬ 
tant implications, both for superconductivity and for pre¬ 
dicted topological order in FeSe/SrTiOs. First, the new 
F band will serve as an essential input for revised FRG 
calculations of the effective low-energy pairing interac¬ 
tion [9]. Second, the modest 140 meV gap we measured 
between filled and empty Gamma bands gives hope that 
inversion of these bands may be achievable, and may lead 
to a predicted topological phase [32]. Finally, our work 
introduces decay length spectroscopy as a general and 
complementary technique to QPI imaging, to map the 
absolute momentum-resolved electronic band structure of 
filled and empty states using STM. We suggest the use of 
these techniques in concert to track the F pocket energies 
in future strain engineering experiments with FeSe. 
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Eigure SI presents the individual energy layers of \g{q^oo)\ that constitute the high-energy portion of Eig. 3(j) of 
the main text. There is one set of rings at -20 meV that grows with increasing energy (branch 1 in Eig. 3(j), red in 
Eig. SI), and another set of higher intensity rings that first appears around 70 meV (branch 3 in Eig. 3(j), blue in 

Fig. SI). 
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FIG. SI. Quasiparticle interference (QPI) imaging, momentum-transfer (q) space, (a-o) Fourier transform amplitudes \g{q,uj)\ 
of conductance maps used in Fig. 3(j) of the main text. All images are displayed with the same color scale range. The red and 
blue arrows denote QPI branches 1 and 3 as discussed in the main text. 

In general, near band edge energies, both quasiparticle interference (QPI) and Bragg (lattice) signals may have 
weight near the reciprocal lattice vectors. Here, we illustrate in Fig. S2 that our data exhibit three characteristic 
distinctions between the Bragg and QPI signals: (1) gr-space extent, (2) intensity, and (3) dispersion. The Bragg 
peaks are only a few pixels wide in the raw data; they have large intensity; and they appear in the same gr-space 
location across all energies. Furthermore, they collapse onto a single pixel after applying drift-correction according 
to the topographic lattice simultaneously acquired (see Refs. [1, 2] for extensive details). On the other hand, the QPI 
signals possess a larger radius in gr-space; they appear at lower intensities; and they disperse/appear/disappear with 
energy. The QPI signal that first appears around 70 meV, corresponding to the emergence of an empty-state band, 
cannot be explained by an abnormally large smearing of the Bragg signal. 

Figure S3 illustrates the azimuthal averaging procedure used to produce Fig. 3(j). Within a large radius of G = 
(27r/a, 0) encompassing the QPI signal, every pixel of \g{q,uj) \ is averaged with those having identical distance values 
Qr to G, then displayed as a function of The resulting plot makes full use of our gr-space pixel resolution, but the 
discrete Qr values are spaced unevenly (Figs. S3(a,b)). We bin each discrete Qr value with its nearest integer pixel 
value to produce the final plot in Fig. S3(c). We additionally note that the first pixel column in Figs. S3(b,c) is nearly 
saturated because the Bragg peak intensity is so much greater than the QPI signal. 

The three bands observed in this work (M electron pocket, F hole pocket, F electron pocket) are quantified by a 
combination of three techniques: (1) QPI |g'(g^,cj)!, (2) decay length spectroscopy dXjduj^ and (3) d‘^IIdV‘^. Fig¬ 
ures S4(a-c) show constant energy cuts of |g(g^,cc;)|, with prominent dispersing peaks due to M electron pocket 
scattering. These dispersing peak positions are fit to a parabola (Fig. S4(d)), which then serves as the guide overlaid 
in Fig. 3(j) of the main text and labeled “branch 1”. Slight deviations from the parabolic fit are present within the 
superconducting gap energy 2A. 
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FIG. S2. Distinct observations of Bragg signals at the reciprocal lattice vectors and QPI signals around the reciprocal lattice 
vectors. (a,b) Conductance maps at two energy layers, -20 meV and 80 meV, along with their Fourier transform amplitudes 
\g{q^u)\ in (c,d). (e,f) After applying a drift-correction algorithm using parameters determined entirely from the simultaneous 
topographic map, the Bragg peaks appear as a non-dispersive, single-pixel entities. (g,h) By examining lower intensities of 
(e,f), QPI signals manifest as dispersive rings with larger qr-space radius. 



FIG. S3. Azimuthal averaging procedure used for visualizing QPI dispersion. (a,b) First, every pixel of \g{q^ijj)\ near G — 
(27r/a, 0) is averaged with those having identical distance values qr to G, then displayed as a function of qr. (c) Second, each 
discrete qr value is binned with its nearest integer pixel value to produce the final plot shown in Fig. 3(j). 


QPI branches 2 and 3 are visible in Figs. S4(a,c) as emerging peaks centered about = 0 that grow in amplitude 
and width away from the Fermi energy Ep. The band edge for branch 2 is difficut to fit, along either the energy or 
qr axes, because its signal overlaps with that of branch 1. However, its hole-like dispersion is evident when one tracks 
the tail end of the dispersing peaks (green pixels in Fig. S3(c) within qr ^ 0.2-0.3 [27r/a]) and observes a change in 























































FIG. S4. (a-c) Azimuthally-averaged line cuts of offset in evenly-spaced energies. Qr is measured relative to G = 

(27r/a, 0). Dispersing peaks from M electron pocket scattering (branch 1 in the main text) are fit to Gaussians (black line 
overlays), and the extracted peak positions are fit to a parabola in (d). The superconducting gap is marked by 2A. 


dispersion direction near -60 meV. (A similar kink in the dispersion is visible at high energies, corresponding to the 
onset of branch 3). We quantify branch 2 by tunneling decay length and dl/dV measurements. Details are discussed in 
the main text, but here in Fig. S5(c,d) the numerical derivatives dXjduj and d‘^IIdV‘^ and their peak fits are explicitly 
shown. We overlay a guide for branch 2 in Fig. 3(j) with its band edge informed by these two measurements (-59 
± 5 meV and -65 ± 3 meV respectively) and its dispersion informed by angle-resolved photoemission spectroscopy 
(ARPES) measurements of the same pocket [3]. 



FIG. S5. (a) Gonstant-^r line cuts of \g{qr,(^)\, offset in evenly-spaced qr values and normalized by their maximum amplitude 
for improved visualization. Dispersing peaks from F electron pocket scattering (branch 3 in the main text) are fit to Gaussians 
(black line overlays with horizontal bars indicating peak positions), and the resulting peak positions are fit to a parabola in 
(b). (c) Energy dependent decay length A(ct;) (red circles), repeated from Fig. 4(b) of the main text, along with its numerical 
derivative (blue circles), (d) dl/dV spectrum (red solid line), repeated from Fig. 4(c) of the main text, along with its numerical 
derivative (blue dashed lines). Gaussian peak fits for band edge extraction are overlaid (black solid lines), and vertical shaded 
bars mark their peak positions. 

To capture the band edge of branch 3, we take constant-g^^ cuts of \g{qri^)\^ and fit the observed peak locations 
to a parabolic dispersion (Figs. S5(a,b)). The resulting fit serves as the guide for branch 3 overlaid in Fig. 3(j). The 
position of this line is further confirmed by the positive energy extrema of dXjduj and d?IldV‘^^ which yield 78 ± 5 
meV and 77 ± 3 meV respectively. 

A summary of the band parameters extracted from these quantitative analyses is given in Table SI. 
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kp [ir/a] 

m*/m 

So [meV] 

eo [meV] 

3(j) 

S4, S5(a,b) 

\9{qr,0j)\ 

0.22 ±0.01 

2.0 ±0.1 

(-60) 

75 ±3 

4(b) 

S5(c) 

dX/dui 

- 

- 

-59 ±5 

78 ±5 

4(c) 

S5(d) 

d^IjdV"^ 

- 

- 

-65 ±3 

77 ±3 


TABLE SI. Summary of band parameters quantified by QPI imaging \g{qr^oj)\^ decay length spectroscopy dXjduj^ and d^IjdV^. 
kp denotes the Fermi wave vector, m* is the effective mass, and so is the band edge. 


Finally, we present additional decay length spectroscopy measurements. Figures S6(a,b) compare X{uj) measure¬ 
ments of single-unit-cell FeSe/SrTiOs taken with a Ptir tip with two different microscopic terminations. While the 
absolute values of A(co’) vary slightly due to differences in tip work function, the steep rises corresponding to the onset 
of F-centered pockets occur at the same energies. Figure S6(c) shows a calibration measurement on polycrystalline 
An exhibiting a flat A(cj), consistent with previous measurements on Au(lll) [4]. 


Tip A - 1 UC FeSe/STO Tip B - 1 UC FeSe/STO Tip C - Poly Au 



__l _I_I_I_ I Q 0.30 _ I_I_I _ ■ Q n 30 h ■ _I_I_I_1_ 

-200 -100 0 100 200 -200 -100 0 100 200 -200 -100 0 100 200 
Energy [meV] Energy [meV] Energy [meV] 


FIG. S6. (a,b) Energy dependent decay length A(ct;) of single-unit-cell (lUC) FeSe/SrTiOa, taken with a PtIr tip with two 
different microscopic terminations (called A and B). (c) A(ct;) acquired on polycrystalline Au. 


* Present address: State Key Laboratory of Low-Dimensional Quantum Physics, Department of Physics, Tsinghua University, 
Beijing 100084, China 
^ jhoffman@physics.harvard.edu 

[1] M. J. Lawler, K. Fujita, J. Lee, A. R. Schmidt, Y. Kohsaka, C. K. Kim, H. Eisaki, S. Uchida, J. C. Davis, J. P. Sethna, 
and E.-A. Kim, Nature 466, 347 (2010). 

[2] 1. Zeljkovic, E. J. Main, T. L. Williams, M. C. Boyer, K. Chatterjee, W. D. Wise, Y. Yin, M. Zech, A. Pivonka, T. Kondo, 
T. Takeuchi, H. Ikuta, J. Wen, Z. Xu, G. D. Gu, E. W. Hudson, and J. E. Hoffman, Nat. Mater. 11, 585 (2012). 

[3] D. Liu, W. Zhang, D. Mou, J. He, Y.-B. Ou, Q.-Y. Wang, Z. Li, L. Wang, L. Zhao, S. He, Y. Peng, X. Liu, C. Chen, L. Yu, 
G. Liu, X. Dong, J. Zhang, C. Chen, Z. Xu, J. Hu, X. Chen, X. Ma, Q. Xue, and X. J. Zhou, Nat. Comm. 3, 931 (2012). 

[4] Y. Zhang, V. W. Brar, F. Wang, C. Girit, Y. Yayon, M. Panlasigui, A. Zettl, and M. F. Crommie, Nat. Phys. 4, 627 (2008). 
















